Evidence for two fluids in the ground state of cuprates 
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We report charge density measurements, using NMR, in the superconducting compound 
(Ca a; Lai_ a ;)(Bai.75_ I Lao.25+a;)Cu30j ; , which has two independent variables x (family) and y (oxy- 
gen). For underdoped samples we find the rate at which holes are introduced into the plane upon 
oxygenation to be family-independent. In contrast, not all carriers contribute to either antiferromag- 
netic or superconducting order parameters. This result is consistent with a two fluid phenomenology 
or intrinsic mesoscopic inhomogeneities in the bulk. We also discuss the impact of weak-chemical- 
disorder on T c . 



Proper counting of holes in the cuprates is essential for 
understanding their properties. Since they have a small 
ratio of coherence length to mean free path, they are con- 
sidered clean superconductors where the carrier density 
n is the same as the superconducting carrier density n s 
in the zero temperature limit 0. In contrast, modern 
experiments in the superconducting part of the phase di- 
agram consistently find strong inhomogeneities in these 
materials 0, 0- Moreover, several theories of cuprate 
superconductivity are based on two fluids comprised of: 
hole-poor and hole-rich regions!^] , bosons and fermions 
H, hot and cooled electrons Q, etc.. This contradic- 
tion leads to the question: do all holes participate in 
the superconducting order parameter? Addressing this 
question requires a simultaneous measurement of n and 
n s . The in-plane 63 Cu(2) nuclear quadrupole resonance 
(NQR) parameter vq is a direct measure of n. n s can 
be extracted from the penetration depth. In this work 
we compare the two numbers for different superconduct- 
ing families. We find that n s is not a universal fraction 
of n, and provide a simple relation between n s and y, 
from parent to overdoped samples. The data support 
two fluids-based theories for the entire phase diagram, 
but constrain the possible division between the fluids. 

Our study is done on the 

(Ca a! Lai_ a ;)(Bai.75-a:Lao.25+x)Cu 3 1 , (CLBLCO) 
compound where each value of a; is a superconducting 
family with its own maximum T c {T™ ax ) ranging 
between 58 and 80 K, as shown in Fig. QJa) @|. By 
choosing a sample with a particular x and oxygen 
level y one can control n and n s independently. This 
allows determination of the rates at which oxygenation 
produces carriers and carriers turn superconducting. 
We evaluate the average vq form 63 Cu(2) nuclear 
magnetic resonance (NMR) measurements. In passing, 
we estimate the chemical disorder from the width, Avq, 
of the vq distribution. n s is obtained from previous 
muon spin rotation (/iSR) measurements 0, • 

The ability of NMR to determine carrier density is 
based on the fact that 63 Cu, with its spin 3/2 nuclei, 
is directly coupled to charge degrees of freedom via the 
electric field gradient (EFG), and vq is a measure of this 
coupling, vq, in turn, depends linearly on the hole den- 



sity |8j according to 
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where A and Vq are doping-independent, but, in princi- 
ple, could be family-dependent. This linear dependence 
was demonstrated for various compounds such as Y123 
0, La 2 - r£ Sr a; Cu04 (La214) 0, and HaBa 2 Cu0 4 +5 
(Hal24) [Hi. Therefore, Eq. [T]and the ability of NMR 
to detect the in-plane copper [Cu(2)] vq selectively will 
allow us to determine the evolution of the in-plane carrier 
concentration and the width of its distribution. 

The measurements were done on powder samples fully 
enriched with 63 Cu. Their preparation is described in 
Ref. 01 • The oxygen content was measured by double 
iodometric titration. The accuracy of this method in the 
enriched CLBLCO is about 0.01. We measured between 
five and seven different samples for each x in the nor- 
mal state at 100 K. The most overdoped sample is a non 
superconducting x = 0.1 compound. The NMR mea- 
surements were done by sweeping the field in a constant 
applied frequency / app = 77.95 MHz, using a 7r/2 - it 
echo sequence. The echo signal was averaged 100, 000 
times and its area evaluated as a function of field. The 
full spectrum of the optimally doped x = 0.4 sample 
(y = 7.156) is shown in the inset of Fig. [2] 

The Cu spin Hamiltonian can be written as |l3l |: 

H/h = -vfl • (1 + K) • H+^[3I2-I 2 +^(I2-I2)], ( 2 ) 

where vi = ( 63 -f/2ir)H, H is a unit vector in the direc- 
tion of the field, K is the shift tensor, and rj is the asym- 
metry parameter of the EFG. In the absence of mag- 
netic field, there is only one transition frequency given 
by / = VQyfl + n/3, so vq cannot be separated from 
n, and the use of the magnetic field is essential. This 
field, applied in the direction 9 and <j) with respect to 
the principal axis of the EFG, lifts this degeneracy, and 
three transition frequencies v m {H,6,Lp) are expected: a 
center line which corresponds to the 1/2 — ► —1/2 tran- 
sition (to = 0), and two satellites which correspond to 
the 3/2 -» 1/2 (m = 1) and -1/2 -> -3/2 (m = -1) 
transitions. Expressions for v m (H, 8, ip) up to second or- 
der perturbation theory in vq for completely asymmetric 
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FIG. 2: NMR spectra of 63 Cu at T = 100 K in optimally 
doped CLBLCO samples with varying x. The inset shows the 
full spectrum of the x = 0.4 compound including contribu- 
tions from Cu(l) and Cu(2). The main figure zooms in on 
the Cu(2) contribution (note the three axis breakers). The 
position of the Cu(2) peaks are shown by dotted lines. 



FIG. 1: (a) The superconducting part of the 
(Ca ;c Lai_ ;c )(Bai.75_ ;c Lao.25+a;)Cu30j ; phase diagram [l|- 
(b) Normalized critical temperatures plotted as a function of 
K(x)Ay, including the glass T g temperature, and corrected 
Neel temperature T^ r (see text) 0. Inset, the normalized 
muon spin relaxation rates taken from Ref. @] as a function 
of oxygen level y for the x = 0.1 and 0.4 CLBLCO samples. 

EFG and shift tensors are given in Ref. 0] . In a powder 
spectrum, where 9 and ip are integrated out, each one of 
the m = — 1 and m = 1 transitions contributes one peak, 
and the m = transition generates two peaks provided 
that rj < 1 [lil Hif. In principle, the bigger vq is, the 
further the peaks are away from each other. Similarly, 
the wider the peaks are, the broader the distribution of 
vq. In CLBLCO as in YBCO, the middle peak is from 
the Cu(l) which has r\ ~ 1. This peak is labeled 1 in the 
inset of Fig. [21 The other four peaks are associated with 
the plane Cu(2) and are labeled 2. 

A zoom on the main features of the Cu(2) signal of 
all four optimally doped samples is depicted in Fig. [21 
(note the three axis breakers). The evolution of the main 
peaks as x increases is highlighted by the dotted lines. It 
is clear that as x decreases the peaks move away from 
each other. This means that vq at optimal doping is 
a decreasing function of x. A more interesting obser- 
vation is the fact that there is no change in the width 
of the peaks, at least not one that can easily be spot- 
ted by the naked eye. This means that the distribution 
of vq is x-independent and that there is no difference 



in the weak-chemical-disorder (WCD) between the opti- 
mally doped samples of the different families. By WCD 
we mean disorder that is not strong enough to "wipe 
out" the contribution of the Cu(2) in its vicinity from 
the spectrum. Thus, WCD is not relevant to the varia- 
tion of T™ ax between the different families. As we shall 
see, this conclusion is supported by more rigorous anal- 
ysis. 

Although wc limit our conclusion to WCD, it is impor- 
tant to place a few limitations on the possible existence 
of strong chemical disorder. The difference in ionic size 
between La and Ba in the Ba layer creates local lattice 
distortions. These impact the apical oxygen 0(4) posi- 
tion, which contributes to v q . However a simple calcu- 
lation based on Eqs. 4 and 5 in Ref. shown that 
the probability of finding an 0(4) out of its ideal place 
increases with x. Therefore, if the 0(4) displacement 
were important for chemical disorder we would have ex- 
pected the x — 0.4 family to have a lower T c than the 
x = 0.1, in contrast to observation. In addition, magic 
angle spinning (MAS) Ca NMR found only one line. This 
ensures the existence of only one Ca site. Also, since Ca 
and Ba have the same valance, and the total amount 
of La is fixed, doping is done only by the chain layer 
oxygen, which is relatively remote from the Cu02 layer. 
Therefore, chemical disorder is expected to be minimal. 
Overall, we cannot find evidence for the existence or im- 
portance of strong chemical disorder. 

The evolution of the main peaks for x — 0.4 as a func- 
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FIG. 3: NMR spectra of 63 Cu(2) at T = 100 K in CLBLCO 
samples with x — 0.4 and varying y. The positions of the 
peaks are shown by dotted lines. 



tion of y is shown in Fig. |3| Here, as y increases the 
peaks moves away from each other, i.e., vq increases as 
a function of doping as expected from Eq.^ Similar data 
for the x — 0.1 family can be found in Ref. [l7| . 

In order to analyze these data more precisely we must 
fit them to a field-swept powder spectrum P(H) which 
is given by 



action of such magnitude in the system. Consequently, 
we added to the numerical evaluation of Eq.|3]a loop over 
200 values of vq drawn from a normal distribution with a 
width Avq. We also added the contribution of the chain 
site, with rj ~ 1. We searched for the best fit to the data 
by x 2 minimization using a simplex code |2(i| . The result 
for the optimal doped samples is shown as the solid line 
in Fig. □ 

The fit is not perfect, mostly because the spectrum 
is not symmetric. Such an asymmetric spectrum is a 
result of correlations between different a priory random 
parameters in the Hamiltonian . For example, if the 
shift tensor K is also non uniform in the sample but its 
values are correlated with the values of vq, the spectrum 
could be non symmetric. Attempts to take this kind of 
effect into account failed due to the enormous increase 
in computing time. We continue the discussion based on 
the best fit we could practically achieve. For a discussion 
on the error bars evaluation see Ref. 0. We could not 
determine, K Zl K x and K y very accurately, and found 
very small r\ for all the samples, ranging from to 0.1, in 
agreement with various estimates for Y123 [j^ . 

The results for fitted vq and Avq are shown in Fig. 0] 
(a) and (b), respectively. From Fig.^Ja) it is clear that 
vq grows linearly with doping in the underdoped side of 
the phase diagram, in agreement with other compounds 
0, Q] , and Eq. ^ It is also clear that Vq in this 
equation is ^-dependent, but this could be attributed to 
NQR base line properties. In contrast, the behavior in 
the overdoped side of the x — 0.1 is surprising since vq 
saturates. As pointed out before in Ref. 17], the added 
holes on the overdoped side no longer go into the planes. 

The most interesting finding is that within the experi- 
mental error, the slope of vq(x, y) in the underdoped side 
is x-independent, as demonstrated by the parallel solid 
lines, and dv Q {x,y)/dy = 7.2(1.4) MHz/Oxygen. This 
means that the rate at which holes are introduced into 
the Cu02 planes, dn/dy, is a constant independent of x 
or y in the underdoped region. Using further the ubiq- 
uitous assumption that the optimal hole density, n opt , at 



f uitous assumption that the optimal hole density, n p , at 

P m (H) oc J M <5[/ app - v m {H, 6, tp)]dtt (3) optimal oxygenation, y opt , is universal, we conclude that 



where M is a matrix element. Usually one can extract 
all the averaged Hamiltonian parameters from the posi- 
tion of the peaks in the spectra 0|- However, we are 
interested in both the parameters and their distribution. 
Therefore, we must account for the entire line shape. For 
this purpose we simulate the line by using a grid in the 
(9, (j>) space, calculate the frequency v m (H, 9, ip) for every 
m, field, and point on the grid, and add 1 to a histogram 
of H when one of the frequencies v m equals / app - The 
matrix elements are taken as unity. This numerical sim- 
ulation approximates P(H) in Eq. To account for the 
peak widths we assume that the main contribution to this 
width, of the order of a few MHz, is from a distribution 
in vq , since the quadrupole interaction is the only inter- 



n(x,y) is a function only of Ay — y — y opt . This re- 
sult was predicted previously [18| based on bond valance 
calculations. 

An implication of this conclusion is that the Presland 
et al. formula 19] T c /T™ ax = 1 - 82.6(n - 0.16) 2 could 
not be correct for CLBLCO. For example, the x = 0.1 
and x = 0.4 samples with T c = have different Ay, and 
according to our experiment different n. However, the 
Presland formula predicts the same n. 

If the conversion from oxygen to holes is universal, but 
different families begin to superconduct at different y val- 
ues (see Fig.nj, it means that the rate at which holes turn 
into superconducting carriers is family-dependent. To 
demonstrate this concept further we focus on the region 
of the CLBLCO phase diagram where T c grows linearly 
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FIG. 4: (a) Nuclear quadrupole resonance frequncy of 63 Cu(2) 
in all the CLBLCO samples extracted from the NMR spe- 
tra as described in the text, (b) The width of the nuclear 
quadrupole resonance frequency distribution of 63 Cu(2) in all 
the CLBLCO samples extracted from the NMR spectra as 
described in the text. 

with doping, as emphasized by the solid lines in Fig.^a). 
This behavior is unique to CLBLCO, which does not 
show a 1/8 dip or plateau in T c . We determine n s from 
/iSR measurement at T — ► 0. The /iSR relaxation rate a 
is a measure of the density of superconducting carriers, 
and a/a opt — n s /n° pt . In the inset of Fig. ^b) we depict 
a/a opt , taken from Ref. 0, as a function of y for two 
extreme families with x = 0.1 and 0.4. In this particular 
region dn s /dy is a constant, which is family-dependent 
and denoted hereafter by K{x). Using the universality 
of dn/dy we find that dn s /dn = K(x)(dn/dy)^ 1 , which 
varies between families. This, again, rules out the possi- 
bility that n s — n in all CLBLCO samples. 

In light of previous work a stronger conclusion could 
in fact be extended to the entire phase diagram of 
CLBLCO. In Fig. Eft)), takcn from Ref - 0> we de P ict 
Tff r ,T g and T c , normalized by T c max of each family. 
Tff r is the Neel temperature after the contribution from 
anisotropies have been divided out. The scaled vari- 



able K(x)Ay, with K = 0.77, 0.67, 0.54, 0.47, for the 
x = 0.1 to 0.4 respectively, collapsed the entire phase di- 
agram into a single function. However, no interpretation 
was given for K(x). The present experiment suggests 
that dn s /dy = K(x) over the entire phase diagram, and 
that n s should be considered as carriers participating in 
both superconducting and antiferromagnetic order pa- 
rameters. 

As for the line width, examination of Fig. EJb) reas- 
sures us of our previous intuition that Auq is not chang- 
ing between the different families, especially for the op- 
timally doped samples. This becomes obvious when con- 
sidering Au Q of YBCOy at T = 100 K, which also has 
v Q = 31 MHz, but Au Q = 0.5 MHz 0. The changes in 
Auq between x = 0.4 and x = 0.1 are minute compared 
to the changes between x = 0.1 and YBCO7. Yet x = 0.4 
has a T™ ax very similar to YBCO7. Thus, if YBCO is 
considered disorder free, the difference in WCD between 
x = 0.4 and x = 0.1 either does not exist or is not rel- 
evant. The same conclusion was reached by MAS Ca 
NMR 0. 

In summary, we find that the rate at which holes are 
doped into the planes when oxygen is added, dn/dy, is 
identical in all families in the superconducting under- 
doped region. In part of this region, the rate at which 
holes contribute to the condensate dn s /dn oc K{x) is a 
family-dependent constant. Based on this and previous 
findings, we conclude that in CLBLCO not all oxygens 
contribute holes to the superconducting and antiferro- 
magnetic order parameters. We also show that weak- 
chemical-disorder, as determined from Auq, is not play- 
ing a role in the T™ ax variations. This reinforces pre- 
vious conclusions that the changes in T™ ax between dif- 
ferent CLBLCO families is caused o nly by variation in 
the in-plane coupling constant J 0, l23| . Finally, since 
u q is not changing in the overdoped side as the oxygen 
level increases, in this side holes are not added to the 
planes. The open question in this study is which physi- 
cal or chemical property of CLBLCO sets K{x). 
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